[X, map] = imread('E:\Seminararbeit\TestImg\storm.png');

A = double(X(:,:,1));


    
    psf = fspecial('motion', mot, deg)
    
    H = conv2(A,psf,'valid')
    
    HD = MyDeconv2DFast(H, psf,size(H,1)-size(psf,1),size(H,2)-size(psf,2),iter)